LA-UR  88-3118 


iiJIC  FitTHV* 


&>■ 


Los  Alamos  National  Laboratory  <s  operated  by  the  University  ot  California  for  the  United  States  Department  of  Energy  under  .ontraci  W-7405  ENG-36 


AD-A201  298 


TITLE  EARLY  TIME  STRUCTURING  AT  VERY  HIGH  ALTITUDES: 

INSTABILITY  MECHANISM,  PROPERTIES  AND  CONSEQUENCES 


AUTHOR(S)  Dan  Winske,  X-l 


submitted  to  Report  prepared  for  the  Defense  Nuclear  Agency 

under  Project  Code  RB,  Task  Code  RC,  Work  Unit  167 

DTIC 

SELECTE 
OCT  1  91988 

% 

H 

By  acceptance  of  this  article  the  publisher  recognizes  that  the  U  S  Government  retains  a  nonexclusive,  royalty-free  license  to  publish  or  reproduce 
the  published  form  of  this  contribution  or  to  allow  others  to  do  so,  tor  U  S  Government  purposes 

The  Los  Alamos  National  Laboratory  requests  that  the  publisher  identify  this  article  as  work  performed  under  the  auspices  of  the  U  S  Department  of  Energy 


FORM  NO  836  «4 
ST  NO  2829  5 '0i 


DISTRIBUTION  STATEMENT  ~A 

Approved  fear  pobttc  ivImhi 
Distribution  Unlimited 


Los  Alamos  National  Laboratory 
Los  Alamos, New  Mexico  87545 

88  10  18  161  - 


Li  i  L.U 


mjagn«afccfci  man 


u.  IMPORT  ilCU«TY  ClAJWOnCK 

UNCLASSIFIED 


AUTNONTY 


AEPO*T  DOCUMENTATION  PAGE 


'*  WTJUCT1VI  MARKING 


ZB.  0»CUM»KAH0W  /  OOWNOBAOwa  SCHtDUU 


4  MRRORMlNG  ORGANIZATION  RfRORT  NUMBIRtSl 

LA'UR-88-3118 


W  KAMI  OS  MRRORMlNG  ORGANIZATION 

Los  Alamos  NaCional  Laboratoi 


Be  AOORfSS  (Ofjr,  km.  an*  OSC9U) 

Los  Alamos,  New  Mexico  87545 


Be.  NAMI  OR  RUNOlNG  '  SPONSORING 
ORGANIZATION 

Defense  Nuclear  Agency 


•C  AOORIIS  COty,  SUM.  1*  9CMJ 

Washington  DC  20305-1000 


I.  MONITORING  ORGANIZATION  RIRORT  NUMItR(S) 


T4.  NAMI  OR  MONITORING  ORGANIZATION 


tv  AOORUStOQr,  Sum.  aa«  Of 


B.  RNOCURBMRNT  tRSTRuRMNT  ©BNTVlCAnON  NUMBIR 


TO  tOuRQ  OR  RURtOMRG  NUMMRS 


■7TT 


task 

NO. 

1  RB 

RC 

ii.  mu  - 

EARLY  TIME  STRUCTURING  AT  VERY  HIGH  ALTITUDES:  INSTABILITY  MECHANISM,  PROPERTIES,  AND 
CONSEQUENCES 


TZ.  OfRSORIAL  AUTHORS) 


D.  Winske 


TlA.  TVBt  OR  RBBOKT 

Technical 


wM9K;1X€YJ:1 


14.  OATl  oo  RUMT  9M. 

88  September  15 


S.  RAG1  COUNT 

52 


IB.  SURRUMCNTAAY  NOTATION 

This  work  was  supported  by  the  Defense  Nuclear  Agency  under  Project  Code  RE,  Task  Code  RC, 
Unit  Work  Code  167,  Unit  Work  Title:  Simulations  and  Modeling  of  HANE/VHANE 


Flute  instability 


Thu  report  consider!  structuring  that  resuui  from  plasms  streaming  at  snb  Aifvenic  speeds  across  an  external 
magnetic  field,  as  might  occur  in  a  eery  high  altitude  nuclear  explosion  (VHANE).  In  previous  reports  we  have 
'  proposed  the  lower  hybrid  drift  instability  enhanced  by  the  deceleration  of  the  plaama  by  tbe  field  produces  tbe  flute 

modes  observed  on  the  surface  of  expanding  laser  produced  plasmas.  We  have  derived  an  appropriate  dispersion 
equation  to  describe  the  properties  of  the  unstable  waves  and  carried  out  particle  simulations  to  show  the  growth  and 
evolution  of  the  instability.  We  review  the  salient  features  of  this  earlier  work  and  then  describe  recent  additions  and 
refinements  to  the  theory  and  simulations.  We  discuss  improvements  to  the  dynamics  of  a  3-D  expansion,  some  initial 
work  on  a  nonlocal  linear  theory,  the  addition  of  finite  electron  beta  effects  to  the  linear  theory,  tbe  scaling  of  the 
wave  properties  with  the  ratio  of  the  ion  gyroradius  to  the  magnetic  confinement  rsdins,  and  a  better  understanding 
of  the  nonlinear  evolution  of  the  instability.  Furthermore,  we  have  considered  nonsymmetric  expansions,  simulations 
of  the  instability  in  3-D,  the  very  weak  limit  of  the  instability,  and  the  effect  of  colliaions.  Next,  the  consequences 
of  these  results  both  for  the  NRL  laser  experiment  as  well  as  other  experiments  and  for  VHANEs  are  considered .  'X' 
comparison  of  the  linear  and  nonlinear  properties  of  the  waves  observed  in  the  simulations  with  those  seen  in  the 
i  User  experiments  is  carried  out  and  issues  which  farther  experiments  might  reeolve  are  ontlined.  Finally,  the  scaling 

_ of  the  results  presented  here  to  the  much  weaker  instability  regime  expected  in  a  VHANE  is  discussed.  _ 
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Abstract 

This  report  considers  structuring  that  results  from  plasma  streaming  at  sub  Alfvenic 
speeds  across  an  external  magnetic  field,  as  might  occur  in  a  very  high  altitude  nuclear 
explosion  (VHANE).  In  previous  reports  we  have  proposed  the  lower  hybrid  drift  instability 
enhanced  by  the  deceleration  of  the  plasma  by  the  field  produces  the  flute  modes  observed 
on  the  surface  of  expanding  laser  produced  plasmas.  We  have  derived  an  appropriate 
dispersion  equation  to  describe  the  properties  of  the  unstable  waves  and  carried  out  particle 
simulations  to  show  the  growth  and  evolution  of  the  instability.  We  review  the  salient 
features  of  this  earlier  work  and  then  describe  recent  additions  and  refinements  to  the 
theory  and  simulations.  We  discuss  improvements  to  the  dynamics  of  a  3-D  expansion, 
some  initial  work  on  a  nonlocal  linear  theory,  the  addition  of  finite  electron  beta  effects  to 
the  linear  theory,  the  scaling  of  the  wave  properties  with  the  ratio  of  the  ion  gyroradius  to 
the  magnetic  confinement  radius,  and  a  better  understanding  of  the  nonlinear  evolution 
of  the  instability.  Furthermore,  we  have  considered  nonsymmetric  expansions,  simulations 
of  the  instability  in  3-D,  the  very  weak  limit  of  the  instability,  and  the  effect  of  collisions. 
Next,  the  consequences  of  these  results  both  for  the  NRL  laser  experiment  as  well  as  other 
experiments  and  for  VHANEs  are  considered.  A  comparison  of  the  linear  and  nonlinear 
properties  of  the  waves  observed  in  the  simulations  with  those  seen  in  the  laser  experiments 
is  carried  out  and  issues  which  further  experiments  might  resolve  are  outlined.  Finally,  the 
scaling  of  the  results  presented  here  to  the  much  weaker  instability  regime  expected  in  a 
VHANE  is  discussed. 


1.  Introduction 

Structure  formation  at  early  times  remains  an  important  issue  for  high  altitude  nuclear 
explosions.  At  low  altitudes  the  coupling  of  the  debris  plasma  to  the  background  occurs 
through  collisional  processes,  and  the  front  edge  of  the  resulting  blast  wave  may  be  subject 
to  Rayleigh-Taylor  type  instabilities  [Brecht  and  Papadopoulos,  1979].  At  higher  (e.g., 
STARFISH)  altitudes,  collisionless  plasma  processes  take  over  with  the  generation  of  a 
high  Mach  number  shock  wave  and  structure  on  the  ion  gyroradius  scale  [Thomas  and 
Brecht,  1986].  At  even  higher  altitudes,  the  so  called  VHANE  regime,  the  Alfven  M»Hb 
numbei  of  the  debris  falls  below  unity  and  only  a  weak  shock,  if  any,  is  produced.  Because 
of  the  lack  of  actual  test  data  in  this  regime,  laboratory  and  space  experiments  as  well 
as  computer  simulations  provide  a  method  to  address  the  question  of  early  time  structure 
formation  at  very  high  altitudes. 

At  STARFISH-like  altitudes  the  debris  gives  up  its  energy  and  momentum  when  it 
has  overrun  an  equivalent  mass  of  air  ions.  At  higher  altitudes,  where  the  air  density  is 
much  less,  the  debris  loses  energy  expanding  against  the  geomagnetic  field.  The  revelant 
distance  over  which  the  debris  comes  to  rest  is  the  magnetic  confinement  radius,  Rb-  As 
the  bomb  plasma  expands,  it  drags  the  magnetic  field  along,  creating  a  diamagnetic  cavity 
of  radius  ~  Rb- 

Diamagnetic  cavities  also  occur  in  space,  both  artifically  and  naturally.  The  artificial 
cavities  result  from  chemical  releases,  such  as  the  AMPTE  (Active  Magnetospheric  Particle 
Tracer  Explorers)  mission  in  which  small  cannisters  of  lithium  and  barium  were  exploded 
in  the  solar  wind,  the  magnetosheath,  and  the  magnetotail  [Krimigis  et  al.,  1982].  Natural 
cavities  are  observed  infrequently  upstream  of  the  Earth’s  bow  shock  [Thomsen  et  al., 

1986] .  They  occur  as  a  result  of  changing  conditions  in  the  solar  wind  that  momentarily 
allow  a  large  fraction  of  the  incoming  ions  to  be  reflected  by  the  bow  shock  and  propagate 
back  upstream  [Thomsen  et  al.,  1988].  Magnetic  field-free  cavities  are  also  formed  around 
the  nucleus  of  a  comet  [Mendis  and  Houpis,  1982],  where  a  large  number  of  emitted  neu¬ 
trals  become  ionized  and  exclude  the  field.  Diamagnetic  cavities  are  also  produced  in  the 
laboratory  by  means  of  laser  produced  plasmas  expanding  against  an  externally  applied 
magnetic  field.  [Okada  et  al.,  1981;  Ripin  et  al.,  1987;  Zakharov  et  al.,  1986]. 

Generally,  “small”  plasmas,  where  small  means  the  ion  gyroradius  pi  based  on  the  ex¬ 
pansion  velocity  and  the  external  magnetic  field  is  comparable  to  or  greater  than  Rb,  tend 
to  show  observable,  field  aligned  instabilities  (structures)  [e.g.,  laser  experiments  (Ripin  et 
al.,  1987;  Okada  et  al.,  1981),  theta-pinch  implosions  (Keilhacker  et  al.,  1974),  plasmoid 
propagation  studies  (Papadopoulos  et  al.,  1988;  Wessel  et  al.,  1988),  AMPTE  releases 
(Bernhardt  et  al.,  1987)],  while  “large”  plasmas  with  Pi  C  Rb  [diamagnetic  cavities  at  the 
bow  shock  [Thomsen  et  al.,  1986]  and  the  cometary  cavity  at  Comet  Halley  [Neubauer, 

1987]  do  not.  We  later  show  that  this  follows  from  the  properties  of  the  instability  involved 
and  thus  conclude  that  structure  of  this  type  is  not  likely  to  be  important  for  the  very 
large  VHANE  plasma. 
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The  experiments  which  do  show  structure  yield  important  information  about  the 
instability  which  generates  it.  In  the  AMPTE  experiments  several  kg  of  barium  were 
released  in  the  Earth’s  magnetotail  and  the  expanding  cloud  was  observed  optically  from 
the  ground  [Bernhardt  et  al.,  1987].  The  plasma  expanded  to  a  final  radius  Rb  ~  p,  ~  200 
km,  while  short  wavelength  (A  ~  40  km)  field  aligned  (flute)  modes  appeared  as  ripples 
on  the  surface  near  the  time  of  maximum  expansion.  The  structures  did  not  appear  to 
move,  suggesting  real  frequency  wr  ~  0,  or  to  change  mode  number  as  the  cloud  began 
to  collapse.  The  surface  waves  observed  in  the  NRL  laser  experiments  [Ripin  et  al.,  1987] 
were  often  more  prominent,  in  that  the  radial  extent  of  the  fingers  was  larger  and  a  smaller 
number  of  modes  was  present.  The  wavelengths  were  still  short,  A  <  pt,  and  again  the 
modes  appeared  to  be  time  stationary.  High  frequency  (~  500  MHz)  noise  corresponding 
to  roughly  the  lower  hybrid  frequency  was  also  detected,  as  has  been  found  at  the  edge  of 
the  barium  cloud  by  the  AMPTE/IRM  satellite  [Gurnett  et  al.,  1986].  Furthermore,  as  the 
external  magnetic  field  was  changed  in  the  laser  experiments,  the  wavelengths  of  the  modes 
did  not  seem  to  be  affected.  The  surface  structures  could  also  been  followed  later  in  time, 
into  their  apparent  nonlinear  stage.  The  tips  of  the  flutes  were  observed  to  freestream 
ahead  of  the  plasma,  tended  to  bend  in  the  electron  cyclotron  sense  in  the  magnetic  field, 
and  sometimes  seemed  to  bifurcate  into  shorter  wavelength  modes.  The  surface  waves  seen 
in  the  Japanese  laser  experiments  [Okada  et  al.,  1981]  were  more  like  the  AMPTE  release 
in  that  they  appeared  at  the  time  of  maximum  expansion  and  then  dissappeared  rather 
quickly.  Although  they  were  also  time  stationary,  unlike  the  NRL  results,  the  observed 
wavenumbers  (k)  had  a  strong  dependence  on  the  applied  magnetic  field,  k  —  Bo  e.  Again, 
lower  hybrid  noise  was  detected.  Finally,  experiments  of  a  similar  nature  have  been  carried 
out  by  Zakharov  et  al.  [1986].  It  is  not  clear  whether  they  observed  the  surface  waves 
visually,  but  for  small  plasmas  (p,  >  Rb)  the  size  of  the  magnetic  cavity  was  significantly 
less  than  Rb,  and  intense  high  frequency  electric  noise  was  detected,  again  suggesting  a 
strong  instability  in  this  regime. 


The  observations  have  stimulated  various  theoretical  investigations  to  understand  the 
instability  that  gives  rise  to  the  structure.  Sydora  et  al.  [1983]  consider  a  diocotron-like 
instability  driven  by  relative  drifts  of  electrons  in  nearby  charge  layers,  although  Akimoto 
et  al.  [1988]  show  that  the  inclusion  of  the  ion  dynamics  in  such  a  model  greatly  reduces 
the  growth  of  the  mode.  The  fact  that  the  instability  appears  even  in  the  absence  of  a 
thin  plasma  ring  is  also  contrary  to  the  assumptions  of  the  diocotron  model.  Peter  et  al. 
[1983]  derive  expressions  for  the  properties  of  a  kinetic  Rayleigh- Taylor  instability,  using 
a  sharp  boundary  analysis.  Hassam  and  Huba  [1987,  1988]  emphasize  the  importance  of 
the  slowing  of  the  plasma  by  the  magnetic  field  and  include  the  ion  inertia  term  in  a  set 
of  modified  MHD  equations.  Characteristics  of  the  instability  in  the  long  wavelength  limit 
appropriate  to  their  model  have  been  derived;  the  non-Rayleigh-Taylor-like  nature  of  the 
mode  has  also  been  stressed.  Okada  et  al.  [1979]  and  later  Winske  [1988]  and  Galvez  et 
al.  [1988]  have  included  both  ion  and  electron  dynamics  and  argued  that  the  instability  in 
question  is  the  lower  hybrid  drift  instability  enhanced  by  the  deceleration  of  the  expanding 
plasma.  The  inclusion  of  electron  effects  give  a  physical,  short  wavelength  cutoff  and  a 
peak  in  the  growth  rate  as  a  well  defined  wavenumber.  Both  Okada  et  al.  [1979]  and  Galvez 
et  al.  [1988]  consider  the  electrostatic  limit  with  finite  temperatures  for  both  the  electrons 
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and  ions.  Winske  [1988]  includes  electromagnetic  corrections  but  takes  the  electrons  as 
cold.  Generally,  the  wavenumbers  at  maximum  growth  calculated  from  these  models  that 
include  electron  effects  are  a  factor  of  3-30  too  large  to  explain  the  observations.  Okada 
et  al.  [1981]  compare  with  their  laser  experiments,  while  Winske  [1988]  and  Galvez  et 
al.  [1988]  with  the  AMPTE  magnetotail  releases.  There  are,  however,  a  number  of  other 
effects  that  can  be  included  in  the  analysis  which  improve  the  agreement.  A  discussion  of 
these  issues  comprises  a  major  portion  of  this  report. 

Furthermore,  the  experiments  and  theory  have  been  supported  by  a  number  of  sim¬ 
ulations  of  plasmas  expanding  against  an  ambient  magnetic  field.  Two-dimensional  elec¬ 
trostatic  particle  simulations  which  show  structure,  but  no  development  of  a  diamagnetic 
cavity,  have  been  carried  out  by  Sydora  et  al.  [1983]  and  Galvez  et  al.  [1988].  Electro¬ 
magnetic  particle  simulations  displaying  cavity  formation  have  been  done  by  Gisler  [1988] 
and  Winske  [1988].  Gisler’s  simulations  were  done  in  the  plane  containing  B  and  show  no 
instability;  Winske’s  runs  were  in  the  plane  perpendicular  to  B  and  show  the  structure 
formation,  although  they  do  not  include  the  dynamics  along  the  field  direction  that  Gisler 
studies.  Other  simulations  have  been  done  with  the  electrons  treated  as  a  fluid:  Sgro  et 
al.  [1988]  including  electron  inertia,  Brecht  and  Thomas  [1988]  with  massless  electrons, 
and  Huba  et  al.  [1987]  with  a  single  fluid  (including  the  Hall  term).  All  three  simulations 
show  the  development  of  the  instability,  even  though  the  geometry  was  different  (Sgro  et 
al.  [1988]  and  Huba  et  al.  [1987]  in  a  slab  geometry  instead  of  the  more  conventional 
cylindrical  expansion)  as  are  the  characteristic  wavelengths  of  the  modes  (determined  to 
some  degree  by  the  cell  size  used  in  the  computations).  The  net  conclusion  is  that  the 
instability  is  a  very  robust  beast.  The  fact  that  it  occurs  in  a  variety  of  simulations  which 
have  different  physical  models  suggests  that  it  can  be  investigated  from  many  theoretical 
points  of  view.  We  stress  that  while  the  various  approaches  tend  to  emphasize  different  as¬ 
pects  and  often  give  a  different  name  to  the  instability,  the  underlying  physical  mechanism 
is  basically  the  same. 

The  purpose  of  this  report  is  to  summarize  DNA  sponsored  research  at  Los  Alamos 
related  to  the  early  time  structuring  problem  at  very  high  altitudes.  We  begin  with  a  review 
of  the  basic  physics  model,  results  from  linear  analysis,  and  the  simulations.  Details  of 
this  work  can  be  found  in  Winske  [1987,  1988]  and  Akimoto  et  al.  [1988].  We  then  discuss 
a  number  of  enhancements  to  the  theory  and  recent  simulation  work  which  bears  on  the 
problem.  Specifically,  we  describe  additions  to  the  basic  model  which  include  an  improved 
treatment  of  the  dynamics  of  the  expansion,  the  addition  of  electron  kinetic  effects  to  the 
linear  theory,  and  some  work  on  a  nonlocal  linear  analysis.  Related  to  the  simulations,  we 
have  done  a  careful  study  of  Pi/Ra  scaling,  an  analysis  of  nonlinear  mode  coalescence,  some 
work  on  asymmetric  expansions,  and  a  preliminary  study  of  additional  physics  that  occurs 
in  three-dimensional  simulations.  We  have  also  considered  the  weak  limit  of  the  instability 
and  included  the  effect  of  collisions  in  both  the  linear  theory  and  the  simulations.  Finally, 
the  application  of  this  work  to  the  laser  experiments  at  NRL  and  VHANEs  is  addressed. 
The  key  question  related  to  the  laser  experiment  involves  the  scaling  of  the  instability  as 
a  function  of  parameters  such  as  the  applied  magnetic  field,  Pi/Ra,  and  the  asymmetry 
of  the  expansion.  We  also  discuss  long  time  effects,  such  as  bifurcation  of  the  tips  of  the 


structures,  and  how  the  properties  of  the  structure  change  with  the  background  pressure 
(i.e.,  collisions).  For  VHANEs  the  principal  issue  involves  the  strength  of  the  instability 
and  its  effect  on  the  size  of  the  diamagnetic  cavity,  the  confinement  of  the  debris  plasma, 
and  the  amount  of  electron  and  ion  heating.  The  results  of  such  studies  can  then  be 
used  as  input  for  longer  time  scale  calculations  that  follow  the  evolution  of  the  disturbed 
atmosphere  for  many  hours  after  the  explosion. 


2.  Review  of  Basic  Instability  Mechanism,  Linear  Theory  and  Simulations 

We  begin  with  a  review  of  the  basic  physics  of  the  instability  and  some  results  of  linear 
analysis  and  simulations.  Much  of  this  work  is  explained  in  more  detail  in  earlier  reports 
[Winske,  1987,  1988].  Later  sections  discuss  recent  refinements  to  this  earlier  work  as  well 
as  its  relation  to  the  structuring  seen  in  the  laser  experiments  and  possible  consequences 
for  VHANEs. 


The  physics  of  the  instability  is  shown  in  Figure  1.  We  consider  a  debris  plasma 
expanding,  either  thermally  or  with  a  directed  velocity,  against  an  ambient  magnetic  field, 
Ba.  For  the  discussion  here,  we  consider  only  the  expansion  across  the  magnetic  field; 
3-D  effects  will  be  discussed  in  Section  3A  and  3G.  On  a  time  scale  short  compared  to 
the  ion  cyclotron  period,  Hit  <  1  (ion  cyclotron  frequency  fl,  =  e50/m,c,  where  e  is  the 
electronic  charge,  m,  is  the  ion  mass,  and  c  is  the  speed  of  light),  the  ions  can  be  considered 
unmagnetized  and  thus  expand  radially  outward.  The  electrons,  tied  to  the  magnetic  field, 
try  to  hold  the  ions  back,  giving  rise  to  a  radially  inward  pointing  electric  field,  Er.  As 
the  ions  on  the  outer  edge  are  slowed  by  Er,  the  ions  inside  catch  up,  compressing  the 
plasma  into  a  thin  shell.  The  electrons  and  magnetic  field  are  also  compressed,  forming  a 
diamagnetic  cavity.  Because  the  electrons  are  magnetized,  they  E  x  B  drift  azimuthally 
relative  to  the  ions.  This  current  (I)  is,  of  course,  in  the  right  sense  to  generate  a  magnetic 
field  opposite  to  the  applied  field  inside.  It  is  the  relative  electron- ion  drift  in  the  shell 
which  is  the  free  energy  source  for  the  instability. 

The  simplest  way  to  carry  out  a  linear  analysis  of  the  perturbations  is  a  local  theory, 
where  the  radial  dependence  of  the  modes  is  neglected.  A  nonlocal  approach  is  described 
in  Section  3B,  where  we  discuss  the  physical  regime  in  which  each  is  valid.  One  considers 
unmagnetized  ions  and  magnetized  electrons  drifting  azimuthally  with  respect  to  the  ions 
with  velocity  VE  in  a  uniform  magnetic  field.  To  get  V#  we  use  the  fact  that  the  radial 
electric  field 

tEr  =  T{Cn  (l) 

where  g  =  —dVir/dt(>  0)  is  the  deceleration  of  the  ions,  7*  is  the  ion  temperature,  and 
c„  =  ~n~1driifdx(>  0)  is  the  density  gradient.  VE  is  thus 


VE  =  - 


cEr 

T 


rrijcg  T\c_ 
eB0  +  eBD€n 


=  Vg  +  V„ 


(2) 
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MAGNETIC 


Figure  1.  Schematic  of  the  structuring  mechanism  showing  the  outward  motion  of  the 
ions,  the  radial  electric  field,  the  electron  E  x  B  drift,  the  resulting  current  (I),  and 
the  diamagnetic  cavity. 


The  second  term  due  to  the  ion  pressure  gradient  is  the  usual  driving  term  for  the  lower 
hybrid  drift  instability  [e.g.,  Davidson  and  Gladd,  1975].  For  most  applications,  where  it 
has  been  studied,  e.g.,  theta  pinches,  the  contribution  of  the  deceleration  is  small  (V9  Kn) 
and  ignored.  As  we  shall  see,  the  situations  of  interest  here,  when  the  instability  generates 
large  scale  structures,  is  characterized  by  Vg  >  Vn.  As  with  earlier  work  [Davidson  and 
Gladd,  1975;  Okada  et  al.,  1979;  Winske,  1988],  we  will  continue  to  use  the  name  lower 
hybrid  drift  instability  when  Vg  ±  0,  in  particular  even  when  Vg  >  Vn. 

The  resulting  dispersion  equation  can  be  written  down  in  various  forms,  depending 
on  what  other  effects  are  included.  For  example,  limiting  the  perturbations  to  being 
electrostatic  and  including  Te  /  0  gives  the  electrostatic  limit  [Eq.  (2)  of  Okada  et  al. 
(1979)  or  Eqs.  (A.1-A.8)  of  Galvez  et  al.  (1988)].  Winske  [1988]  (Eq.  20)  includes 
electromagnetic  effects  but  Tt  =  T{  =  0.  This  latter  form  is  convenient  for  making  analytic 
approximations  and  in  showing  that  the  instability  persists  in  the  long  wavelength  (k  — ►  0) 
limit  where  charge  neutrality  and  massless  electrons  are  assumed.  The  scaling  that  results 
in  this  case  is  equivalent  to  that  of  Hassam  and  Huba  [1987]  and  explains  why  the  instability 
occurs  in  hybrid  simulations  with  mt  =  0  [Brecht  and  Thomas,  1988]  as  well  as  in  single 
fluid  Hall  MHD  codes  [Huba  et  al.,  1987]  with  maximum  growth  usually  occurring  at  the 
shortest  allowable  wavelengths  in  the  calculations. 

Figure  2  displays  solutions  of  the  linear  dispersion  equation  [Eq.  (28)  of  Winske  (1988)] 
for  Vg  =  0,  corresponding  to  the  usual  lower  hybrid  drift  instability,  for  Te  =  0.  (Other 
parameters  are:  enc/ui  =  10,  ui 2  =  47rn,e2/mi,  /?,  =  87rn,Ti/B2  =  0.2  so  that  Vn/vA  =  1, 
=  jB2/47rnim,  with  mt/me  =  1836  and  vA/c  —  1/200).  The  real  (wr,  solid  curve)  and 
imaginary  (7,  dashed  curve)  parts  of  the  frequency  are  plotted  versus  wavenumber.  For 
convenience,  the  frequencies  are  normalized  in  terms  of  D,  and  k  in  terms  of  kc/u;,.  Growth 
of  the  instability  persists  down  to  k  — ♦  0,  where  7  ur,  but  peaks  at  shorter  wavelengths, 
kc/w,  ■ '  150,  implying  kc/ue  ~  4,  with  7  ~  wr  ~  30R*  ~  0.8u >lh,  where  the  lower  hybrid 
frequency  is  defined  as  =  w2/(l  +w2 /H2)  ~  nen,-.  In  contrast,  Figure  3  shows  results 
(with  same  scales)  for  the  same  parameters  except  that  Vg/v *  =  3.  Maximum  growth  of 
the  instability  is  now  larger  (about  a  factor  of  two)  and  occurs  at  smaller  wavenumbers 
(about  a  factor  of  two),  with  again  wr  ~  7  at  maximum  growth  and  wr  <  7  at  smaller  k. 

Analytic  expressions  for  7  and  k  corresponding  to  maximum  growth  can  also  be  derived 
[Winske,  1987,  1988],  For  example,  assuming  Vg  »  Vn  one  finds 


*  =  ULH{tn/g)l/2 


and 

QJ  _  (tnC/Wj)1/2  1  +  ty/3 
WIH  “  {Vg/vAy/«  2 


Relating  g  to  initial  conditions  assuming  a  2-D  expansion  [Winske,  1988]  yields 


1/  c  n* 
9  =  VD - 


(3) 

(4) 

(5) 
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r  hybrid  drift  instability  (’ V9  =  0, 
:)  and  imaginary  ('),  dashed  curve) 


but  Vj/vyi  =  3,  showing  enhancement  of  the  instability  when 
d  by  the  magnetic  field. 


Figure  3.  Same  as  Figure  2, 
the  plasma  is  decelerate' 


where  Vp  is  the  initial  expansion  velocity,  r0  is  the  initial  plasma  radius,  and  u>c  is  the  ion 
plasma  frequency  using  the  initial  density.  Then  Eq.  (3)  becomes 


fnroc\l/2/^g\l/2/m*\I/4 

VD  1  W 


an  expression  that  will  be  checked  by  simulations,  shortly. 


(6) 


Next,  some  results  of  simulations  that  confirm  this  basic  picture  will  be  shown.  Fig¬ 
ure  4  displays  snapshots  of  ion  density  contours  from  a  two-dimensional  electromagnetic 
particle  simulation  at  various  times,  similar  to  those  described  in  Winske  [1987,  1988]. 
The  plasma  initially  forms  a  dense  column  (r0  =  2 c/ue)  which  expands  with  a  velocity 
Vp  =  0.005c  into  a  vacuum.  The  ambient  magnetic  field  (magnitude  given  by  fle/we  =  0.4) 
is  perpendicular  to  the  plane  of  the  simulation.  The  ions  have  mass  m,/me  =  1600  and 
thermal  speed  =  (2 Tj/m,)1/2  =  0.005c;  the  electron  thermal  speed  is  ve  =  0.05c.  As 
the  plasma  expands,  the  deceleration  tends  to  compress  the  cloud  into  a  shell  on  which 
flutes  appear.  By  WL/rt  =  4,  the  expansion  stops.  At  later  times,  most  of  the  ions  start 
to  recompress,  although  the  plasma  on  the  outer  edge  continues  to  move  radially  outward. 
By  counting  ripples  on  the  surface,  one  finds  that  the  short  wavelength  modes  appeared 
to  have  coalesced  into  longer  wavelength  structures  at  later  times.  Although  there  is  not 
much  azimuthal  motion  of  the  perturbations,  we  later  show  that  the  frequency  of  the  short 
wavelength  modes  is  ~  ulh- 

While  the  results  look  qualitatively  similar  to  runs  presented  in  the  earlier  reports, 
there  are  some  quantitative  differences.  Before,  the  modes  tended  to  be  short  in  wave¬ 
length,  at  about  the  shortest  resolvable  by  the  grid.  At  later  times  there  was  more  diffusion 
across  the  field  and  little  evidence  for  coalescence.  The  present  runs  generally  use  hotter 
plasmas,  which  tend  to  increase  the  wavelength  and  slow  the  growth  of  the  instability.  As 
a  consequence,  the  waves  of  interest  are  well  resolved  and  the  plasma  reaches  maximum 
expansion  radius  without  strong  diffusion  across  the  magnetic  field.  Recompression  and 
coalescence  can  then  proceed.  Earlier  work  [Winske,  1988]  also  included  a  background 
plasma;  its  absence  makes  little  difference  on  the  development  of  the  instability. 

The  scaling  of  the  wavelengths  of  the  modes  observed  early  (“linear”)  and  later  (“non¬ 
linear”)  in  the  simulations  can  be  checked  with  linear  theory,  as  a  function  of  the  magnetic 
field  strength  and  the  ion  to  electron  mass  ratio.  Figure  5  compares  wavenumbers  in  terms 
of  kc/ue  normalized  to  RB,  versus  the  applied  magnetic  field  (fie/u>e)  [left  panel]  and 
m,/me  [right  panel].  The  squares  correspond  to  measurements  of  k  when  the  instability 
first  appears,  the  circles  to  much  later  in  the  runs.  The  solid  curves  correspond  to  Eq. 
(6);  good  agreement  with  both  Cie/u>e  and  m</me  scaling  predicted  by  the  equation  is 
seen.  Nonlinearly,  there  is  very  weak  (if  any)  dependence  on  He/we  (as  seen  in  the  NRL 
experiments)  and  m,/me. 

As  mentioned  earlier,  the  instability  has  been  observed  in  other  kinds  of  simulations 
that  have  been  carried  out  at  LANL  and  elesewhere.  Galvez  et  al.  [1988]  have  only 
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^LHt  =  0 


ujlhI  =  4 


WLH  t  —  2 


VLHt  =  ^ 


^LHt  =  3 


ULHt  =  8 


Figure  4.  Results  of  2-D  electromagnetic  particle  simulations;  ion  density  contours  at 
various  times  showing  shell  formation  and  instability  growth  and  evolution. 


Figure  7.  Results  of  2-D  hybrid  simulations  showing  ions  in  space  at  various  times  [from 
Sgro  et  al.,  1988]. 


electrostatic  perturbations  (i.e.,  no  diamagnetic  cavity  forms),  start  with  a  large  plasma 
r0  >■  pi  and  allow  only  a  thermal  expansion.  In  this  case  Vg  <  Vn,  but  the  instability 
develops  nevertheless.  An  example  of  such  a  calculation  is  shown  in  Figure  6,  where  the 
ions  and  electrons  in  configuration  space  are  shown  at  various  times.  In  constrast,  Figure 
7  displays  the  ions  in  x-y  space  for  a  plasma  slab  expanding  against  a  magnetic  field  [Sgro 
et  al.,  1988].  The  calculations  were  carried  out  using  a  hybrid  (fluid  electrons  with  me  ^  0, 
particle  ions),  again  stressing  that  the  instability  persists  without  kinetic  electron  effects 
and  that  it  does  not  depend  on  the  plasma  having  a  cylindrical  shape. 


3.  Recent  Enhancements  to  the  Theory 

This  section  treats  a  number  of  improvements  to  the  basic  theory  and  further  simu¬ 
lation  studies  of  early  time  structuring.  The  discussion  is  divided  into  nine  subsections: 
(A)  dynamics  of  the  plasma  expansion,  (B)  progess  on  a  nonlocal  theory,  (C)  linear  theory 
with  finite  /?e,  (D)  Pi/RB  scaling,  (E)  mode  coalescence,  (F)  asymmetric  expansions,  (G) 
3-D  simulations,  (H)  weak  instability  limit,  and  (I)  collisions  with  the  background  plasma. 
We  treat  each  of  these  topics  in  turn,  with  further  discussion  of  their  implications  for  the 
observations  in  Section  4. 


A.  Dynamics  of  Plasma  Expansion 

The  standard  approach  to  determine  the  maximum  cavity  size  and  deceleration  is  to 
consider  the  pressure  balance  relation: 

\nm,Vl  =  i/fmy  +  |£  V  (7) 

where  N  is  the  number  of  debris  ions,  VB  is  their  expansion  velocity,  t;  =  dr/dt,  and 
where  in  2-D  V  =  irr2L  while  V  =  47rr3/3  in  3-D,  and  take  a  time  derivative  to  obtain 
g  =  —dv/dt.  At  maximum  expansion  (v  =  0) 


V2 

Smax  =  a—  (8) 

where  a  =  1  for  2-D  and  a  =  3/2  for  3-D,  while  RB  is  obtained  by  evaluating  (7)  at  r  =  Rb. 
Expressing  in  terms  of  Vg  =  g/Ui  and  =  Alfven  speed  assuming  the  debris  ions  are 
spread  to  uniform  density  [TV  =  nBV{RB)\ ,  i.e.,  v\  =  =  Vj  at  r  =  RB ,  one 

finds 

Vg/VA  =  aPi/RB  (9) 

While  the  2-D  calculations,  which  assume  the  magnetic  field  is  always  transverse  to  the 
expansion,  makes  sense  physically,  the  corresponding  3-D  version  is  obviously  a  simplified 
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extrapolation  that  is  valid  only  to  some  degree.  More  correct  is  to  assume  that  the  plasma 
expands  along  the  magnetic  field  unimpeded  by  B,  as  simulations  in  which  motion  along 
B  is  included  show  [Gisler,  1988].  One  can  obtain  a  time  dependent  version  of  (7)  for  a 
better  description  of  the  dynamics,  taking  V  =  irr2(t)L(t)  with  L(t )  =  L0  +  V[)t  [Gisler  and 
Lemons,  1988].  In  this  case  the  solution  is  more  complicated  and  involves  Airy  functions. 
One  finds  from  the  solutions  that  the  time  of  maximum  deceleration  is  slightly  greater 
than  the  time  of  maximum  extent  (/2m)  ,  with 


{YqI^a)  max  —  &Pi/Rb 

(10) 

with  a  ~  2.3  and 

with  Rb  defined  using  (7). 

{Rm  /  Rb)iiizx  0.88 

(11) 

Thus,  the  effective  deceleration  is  some  50%  larger  than  expected  from  the  more  simple 
minded  3-D  model.  Using  a  larger  g  in  the  linear  theory,  of  course,  gives  maximum  growth 
rates  that  are  larger  and  occur  at  somewhat  longer  wavelength  (cf.  Figure  3). 


B.  Nonlocal  Theory 

We  have  also  carried  out  a  nonlocal  analysis  of  the  lower  hybrid  drift  instability  in 
order  to  compare  with  the  local  theory.  The  calculations  are  done  in  the  electrostatic  limit 
with  cold  electrons  and  ions.  As  before,  the  electrons  are  subject  to  an  E  x  B  drift,  Ve, 
on  the  surface.  Now,  however,  we  make  no  assumptions  about  the  frequency  u  and  we 
assume  a  sharp  boundary  density  profile:  n(x)  =  n„,  x  <  5,  n(x)  =  0,  x  >  6,  while  the 
electrostatic  potential  6<f>  ~  6<j>(x)ctkv~wt .  The  details  are  found  in  Lemons  [in  preparation, 
1988];  here  we  merely  summarize  the  major  results  and  conclusions. 

With  this  model  one  can,  as  before,  carry  out  a  local  analysis  at  6<t>{b),  obtaining  the 
dispersion  equation 


1  + 


+ 


e«w 


n2(i-u>vn?)  knew(i  -  w2/nf) 


(12) 


where  u  —  w  —  kVs.  If  one  assumes  w/fle  <  1,  one  obtains  the  usual  lower  hybrid  drift 
instability  dispersion  equation  in  the  electrostatic  (fluid)  limit. 

To  do  a  nonlocal  theory  requires  matching  solutions  for  6<t>(x)  for  x  <  6  and  x  >  b  at 
x  =  6.  This  is  done  by  integrating  Poisson’s  equation  from  b  —  e  to  b  +  c  and  letting  e  — ♦  0, 
assuming  6<f>  is  continuous,  and  obtaining 


2  + 


n?(i 


+ 


neu>(i  -  u>2/n D 


-  ^  =  o 


(13) 
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At  this  point  one  needs  to  be  careful  not  to  throw  out  the  w2/fl2  terms  right  away,  but 
instead  combine  the  second  and  third  terms,  obtaining 


2  + 


w: 


n«w(i  -  w/ne) 


-  ^  =  ° 


(14) 


In  this  form  Eq.  (14)  looks  like  the  local  equation  (12)  except  that  1  +  w2/ 02  — ►  2  and 
en/k  — »  1  in  the  third  term.  (If  one  throws  out  the  w2/fl2  terms  in  (13)  first,  the  first 
term  becomes  2  +  w2/fl2.)  In  the  long  wavelength  ( k  — ►  0)  limit  these  differences  do  not 
matter,  as  one  neglects  these  terms  regardless  and  solves  the  remaining  quadratic  equation 
to  obtain 

w  *  y  +  i(kvEni)1/2  (is) 

in  agreement  with  earlier  work  [Peter  et  al.,  1983],  but  contrary  to  Hassam  and  Huba  [1987], 
In  this  form  the  growth  rate  looks  like  that  of  a  Rayleigh-Taylor  instability.  7  ~  fc1/2.  It 
is  also  interesting  to  note  that  in  the  high  frequency  limit  u>/fie  >  1,  one  does  not  obtain 
maximum  growth  due  to  the  lower  hybrid  drift  instability  from  (14),  rather  the  much 
shorter  wavelength  Buneman  instability  dominates,  u>  ~  we(me/2m,)1/3(l  +  i\/3)/2.  On 
the  other  hand,  if  one  throws  out  the  w2/fi2  terms  in  (13)  from  the  start,  one  obtains  a 
dispersion  equation  which  looks  like  the  local  theory  (12)  [Batchelor  and  Davidson,  1976]. 
Numerical  solutions  obtained  by  Batchelor  and  Davidson  [1976],  which  give  the  lower 
hybrid  drift  instability,  however,  are  inconsistent  in  that  u>/fie  is  not  small. 

The  conclusion  from  this  analysis  is  that  one  has  to  be  very  careful  in  applying  either 
the  local  or  nonlocal  theory  to  the  problem.  For  short  wavelength  modes,  where  it  ~  cn, 
the  local  theory  is  most  appropriate.  The  nonlocal  analysis  can  give  inconsistent  results:  it 
either  gives  maximum  growth  at  frequencies  much  higher  than  the  lower  hybrid  frequency, 
|  yielding  solutions  at  wavelengths  much  shorter  than  the  density  scale  length  or  gives  roots 

that  violate  w/fle  <  1.  Such  short  wavelength  modes  are  quite  sensitive  to  the  very  unre¬ 
alistic  sharp  boundary  that  has  been  assumed.  On  the  other  hand,  for  the  long  wavelength 
modes  the  sharp  boundary  approximation  is  valid,  as  the  broad  radial  eigenfunctions  expe¬ 
rience  the  boundary  only  in  an  average  sense,  while  the  local  approximation  in  this  regime 
I  is  likely  to  be  less  accurate.  Thus,  in  this  limit  Eq.  (15)  is  a  better  approximation  to  the 

growth  rate  than  is  the  7  ~  k  scaling  of  Hassam  and  Huba  [1987]. 


C.  Finite  /?e  Effects 

The  linear  analysis  in  Winske  [1988]  has  been  improved  to  include  finite  electron 
temperature  effects.  The  addition  of  such  effects  is  a  straightforward,  but  tedious  procedure 
[Davidson  et  al.,  1977,  Zhou  et  al.,  1983].  The  method  followed  here  is  that  of  Zhou  et 
al.  [1983],  where  the  dispersion  equation  is  derived  in  some  detail.  An  example  of  how 
such  effects  change  the  linear  properties  is  found  in  Figure  8,  where  the  growth  rate  is 
plotted  versus  wavenumber  for  two  values  of  Tt.  The  parameters  are  the  same  as  in 
Figure  3;  the  curves  are  drawn  on  the  same  scale  for  easy  comparison.  At  low  values 
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Figure  8.  Growth  rates  versus  wavenumber  for  the  lower  hybrid  drift  instability  with 
deceleration  including  finite  0e  effects  (Parameters  same  as  Figure  3). 


of  Te  ( 0e  =  0.2, Te  =  T{  in  the  figure),  the  inclusion  of  finite  electron  temperature  has 
three  effects.  First,  because  Ve  is  increased,  the  maximum  growth  rate  is  slightly  larger. 
Second,  because  of  finite  Te,  k  corresponding  to  maximum  growth  increases  slightly.  Third, 
at  small  k,  there  is  an  enhancement  in  the  growth  rate  for  Vg  ^  0  at  modest  /3e.  At  higher 
/?e,  this  enhancement  at  small  k  disappears,  the  wavenumber  at  maximum  growth  increases 
slightly,  but  the  overall  growth  rates  are  reduced. 

The  fact  that  k  increases  with  (3e  may  seem  contrary  to  Figure  21  of  Winske  [1988], 
which  shows  k  at  maximum  growth  decreasing  with  /?.  There,  however,  T-  w as  also  charg¬ 
ing  to  keep  Ti  =  Te,  and  the  decrease  in  k  was  primarily  due  to  an  increase  in  T,.  However, 
for  sufficiently  large  (3e  the  short  wavelength  modes  are  suppressed  [Drake  et  al.,  1983], 
as  also  shown  in  Figure  21  of  Winske  [1988] .  This  effect  is  due  to  strong  damping  of  the 
electrons  by  the  magnetic  field  gradient.  One  can  easily  show  from  Ampere’s  law  that 
VB/B  ~  /?y2(Vf;/t;e)we/c.  The  bottom  line  is  that  finite  /5e  effects  play  a  minor  role  and 
can  be  neglected  in  favor  of  the  simpler  Te  =  0  dispersion  equation,  unless  /?e  »  1. 


D.  Pi/Rs  Scaling 

As  mentioned  in  the  Introduction,  various  experiments  suggest  that  Pi/Rs  —  i°n 
gyroradius/magnetic  confinement  radius  is  an  important  parameter.  The  structures  that 
appear  on  the  surface  of  expanding  plasmas  are  larger  and  more  prominent  for  pi/RB  >  1, 
as  occur  in  the  NRL  laser  experiment  [Ripin  et  al.,  1987].  Zakharov  et  al.  [1986]  suggest 
that  strong  cross-field  diffusion  occurs  in  this  limit,  as  discussed  later  in  this  subsection. 
Furthermore,  plasmoid  propagation  studies  show  that  Px/Rb  ~  1  (where  Rb  is  now  the 
plasma  radius)  separates  different  physics  regimes  [Papadopoulos  et  al.,  1988;  Wessel  et  al., 
1988].  Again,  for  Pi/Ra  >  1»  i.e.,  when  the  ions  are  essentially  unmagnetized,  anomalous 
effects  allow  the  plasmoid  to  propagate  more  effectively  across  the  magnetic  field. 

Figure  9  shows  the  results  of  a  number  of  2-D  electromagnetic  particle  simulations 
where  Px/Rb  was  varied,  by  changing  m,/me  and  ne/we,  keeping  the  magnetic  confinement 
radius  Rb  ~  (mt/m,.)1/2u;e/ne  =  constant  (in  2-D)  [note  p,  ~  (m,/me)we/ fle].  In  each 
case  the  simulations  are  shown  at  the  same  instant  of  time  u>int  =  8;  contours  of  the 
ion  density  are  plotted.  A  number  of  features  are  clearly  evident.  For  large  Px/Rb  the 
plasma  compresses  into  a  thinner  shell,  the  instability  appears  to  develop  quicker,  and  the 
modes  have  longer  wavelength.  In  the  first  two  cases,  indeed  the  instability  develops  so 
rapidly  that  one  sees  the  nonlinear,  coalesced  structures.  And  in  the  smallest  Px/Rb  case, 
the  plasma  has  not  yet  reached  its  maximum  expansion.  In  addition,  for  larger  p%/Rb  the 
final  size  of  the  plasma  is  smaller,  the  ions  are  better  confined,  recompression  of  the  central 
core  is  more  evident,  and  coalescence  to  longer  wavelengths  is  more  likely.  A  number  of 
these  effects  can  be  inferred  from  linear  theory.  For  example,  because  pi/ Rb  ~  Vg/vA,  the 
deceleration  is  more  rapid  and  as  a  consequence,  compression  into  a  thinner  shell  results. 
In  addition,  larger  VgjvA  implies  a  larger  growth  rate  and  longer  wavelengths  (Figure  3). 
Figure  10  reconfirms  this,  by  comparing  the  observed  k’s  in  the  linear  (early  time)  and 
nonlinear  (later  time)  regimes  as  a  function  of  Px/Rb  with  linear  theory.  Overall,  there  is 
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Pi/RB  =  7.1 


Pi/ Rb  =  0.8 


Pi/Re  =  1.7  =  0.2 


Figure  9.  Results  of  2-D  particle  simulations  showing  ion  contour  plots  at  one  time  (u >uit  = 
8)  for  several  runs  in  which  Rb  =constant,  but  pi/Rs  is  varied. 


Figure  10.  Results  from  linear  theory  (solid  line)  showing  wavenumber  at  maximum  growth 
versus  Pi/RB\  closed  circles  (open  circles)  correspond  to  the  simulation  results  at  early 
(late)  times. 


a  good  match;  the  nonlinear  k’s  suggest  stronger  coalescence  at  larger  Pi/Rp-  While  the 
stronger  (i.e.,  larger  free  energy)  expected  at  larger  Pi/ Rb  also  explains  the  larger  size  of 
the  structures,  the  fact  that  the  final  radius  of  the  bulk  of  the  plasma  is  smaller  is  less 
intuitively  obvious. 

The  results,  however,  are  consistent  with  the  experiments  of  Zakharov  et  al.  [1986], 
who  argue  as  follows.  The  plasma  radius  as  a  function  of  time  consists  of  two  terms,  one 
due  to  expansion,  the  other  due  to  diffusion  of  the  plasma  through  the  field: 


r(i)  =  Vpt  -  Ct 2 


(16) 


The  diffusion  is  so  strong  that  the  final  radius  Rm  <  Rb  so  that  the  expansion  goes  just 
as  Vpt  rather  than  as  Rb  sin(Vpt / Rb)  [Winske  (1988)].  Here,  Ct 2  =  (Dt)1/2,  where  D 
is  the  diffusion  coefficient  ( D  =  uc2 /oj2)  and  v  ~  anomalous  collision  frequency  = 
with  f  =  constant.  The  plasma  reaches  its  maximum  radius  Rm  when  dr/dt  —  0,  i.e.,  at 
Im  —  Vpj2C  and 

Rm  =  V-k  =  TT-^r—r,  (17) 


4C  4  (spxIRb)1'2 


Figure  11  plots  Rm/Rb,  a s  measured  in  the  simulations,  versus  Pi/ Rb-  Qualitatively 
consistent  with  Zakharov  et  al.  [1986],  when  pi/Rs  <  1,  Rm/Rb  ~  1  (in  practice  this  is 
used  to  more  accurately  determine  RB  when  thermal  effects  are  included).  On  the  other 
hand,  when  Pi/Rp  >  1,  Rm  <  Rb',  a  straight  line  through  the  points  gives  Rm/Rb  ~ 
(pi/ Rb)~0'3,  which  is  a  weaker  dependence  than  Eq.  (17).  The  derivation  of  (16),  however, 
assumes  a  3-D  expansion  so  that  N  ~  nr(t)3  ~  n(t)(Vpt)3  and  a  constant  resistivity.  In 
fact,  the  simulations  are  2-D,  which  implies  N  ~  f2.  In  this  case  Eq.  (17)  is  slightly  less 
elegant,  because  the  last  term  ~  t3/2;  the  same  method  of  solution  now  yields  Rm/Rb  ~ 
{Pi/Rs)-1-  Better  agreement  occurs  assuming  that  N  ~  t2  but  f  ~  t2  (i.e.,  the  waves  grow 
in  time  and  hence  the  anomalous  collision  frequency  increases  in  time);  then  Rm/Rb  ~ 
(pi/ Rb)~1'Z'  Alternatively,  one  can  argue  that,  consistent  with  the  usual  lower  hybrid 
drift  instability  [Winske  and  Liewer,  1978],  u  is  not  constant  but  varies  as  (pi/ Rb)~1^z, 
then  (17)  yields  Rm/Rb  ~  (/>»'/ Rb)~1^3 •  None  of  these  fixes  is  entirely  satisfactory. 
However,  the  basic  idea  that  when  the  instability  is  strong,  diffusion  of  the  plasma  by  the 
waves  eventually  limits  the  cavity  size  (Rm)  to  be  less  than  the  classical  limit  (Rb)  seems 
physically  reasonable  and  supported  by  the  simulations. 

Two  other  additional  points  should  be  made.  First,  one  can  estimate  when  the  dif¬ 
fusion  picture  should  break  down.  The  usual  plasma  expansion  model  (Sec.  3A)  shows 
that 

r(t)  ~  Rb  sm(VDt/PB)  ~  Vpt  -  (VDt)3/3R2B  (18) 

We  have  ignored  the  second  term  in  favor  of  a  diffusion  term  ~  Ct2.  Obviously,  if  the 
collision  frequency  is  too  small,  the  f3  dynamic  correction  term  in  (18)  determines  the 
radius,  not  the  diffusion.  Again,  one  can  show  that  diffusion  dominates  for  p,  >  Rb- 


Second,  even  though  the  magnetic  field  cavity  grows  to  maximal  size,  the  outer  edges 
of  the  plasma  continue  to  expand  outward.  This  is  already  somewhat  evident  from  the 
simulation  shown  in  Figure  4.  To  see  the  magnetic  cavity  for  this  case  in  more  detail  and 
its  relation  to  the  position  of  the  plasma,  Figure  12  displays  cuts  in  x  through  the  center  of 
the  simulation  region,  showing  the  ion  density  and  the  magnetic  field  Bz  at  various  times. 
One  sees  that  a  well  defined  cavity  forms  at  early  times,  the  plasma  compresses  into  a  thin 
shell,  and  the  cavity  and  plasma  are  of  comparable  size.  At  later  time,  the  cavity  starts 
to  collapse  along  with  the  inner  portion  of  the  plasma,  while  the  other  edge  of  the  plasma 
remains  moving  outward.  This  corresponds  to  the  freestreaming  of  the  flutes  seen  in  the 
NRL  experiment  [Ripin  et  al.,  1987] . 

E.  Mode  Coalescence  and  Nonlinear  Effects 

The  simulation  shown  in  Figure  4  as  well  as  the  larger  pj/Rs  cases  in  Figure  9  indicate 
that  longer  wavelength  modes  prevail  at  late  times.  Several  explanations  for  this  effect  are 
possible.  First,  because  the  longer  wavelength  modes  have  smaller  linear  growth  rates, 
they  will  take  longer  to  grow  to  large  amplitude.  Second,  because  the  longer  wavelength 
modes  appear  during  or  slightly  after  maximum  expansion,  when  the  deceleration  is  larger, 
the  longer  wavelengths  may  be  a  linear  consequence  of  a  larger  effective  g  (cf.  Figures  3 
and  10).  Third,  the  longer  wavelength  modes  may  result  from  nonlinear  coupling  of  the 
shorter  wavelength  modes.  It  is  known,  for  example,  that  the  wavelengths  associated  with 
the  lower  hybrid  drift  instability  in  the  post  implosion  phase  of  a  theta  pinch  are  typically 
about  a  factor  of  two  longer  than  those  calculated  from  linear  theory  at  maximum  growth 
[Fahrbach  et  al.,  1981].  Drake  et  al.  [1984]  have  constructed  a  mode  coupling  theory  whose 
numerical  solutions  give  quantitative  agreement  with  those  observations.  The  theory  is 
based  on  the  idea  that  when  the  waves  grow  to  sufficient  amplitude,  the  nonlinearity  allows 
the  frequency  mismatch  to  be  overcome  and  coupling  of  the  pump  wave  at  wavenumber 
kQ  to  daughter  waves  at  k  ~  kQj 2  can  proceed.  Although  the  theory  has  been  constructed 
in  the  weak  drift  regime  with  Vg  =  0,  which  is  not  directly  applicable  to  our  situation,  the 
physics  of  the  decay  may  persist  anyway. 

Figure  13  displays  the  results  of  a  number  of  simulations  run  with  various  parame¬ 
ters  showing  the  amount  of  coalescence,  as  determined  by  the  ratio  of  the  mode  number 
observed  during  the  linear  growth  phase  (mx,)  to  that  observed  in  the  later  phase  (ttisl), 
versus  Pi/Ra •  Several  points  can  be  inferred  from  the  plot.  First,  the  amount  of  coales¬ 
cence  for  all  of  the  runs  is  between  two  and  three,  with  an  average  value  of  2.5.  Second, 
there  is  no  strong  variation  with  Pi/Rb  •  Third,  there  is  considerable  scatter  in  the  data, 
which  in  part  reflects  the  difficulty  in  determining  mi  and  mxvL,  and  also  the  fact  that 
the  parameters  vary  over  a  wide  range.  The  wide  scatter  may  obscure  any  weak  scaling, 
and  may  also  suggest  that  Pi/Ra  is  not  the  best  parameter  to  characterize  the  coalescence. 
The  overall  results,  however,  seem  in  qualitative  agreement  with  the  Drake  et  al.  [1984] 
theory. 


Figure  12.  Cuts  through  x  of  the  simulation  of  Figure  4,  showing  the  ion  density  profile 
and  the  diamagnetic  cavity  at  various  times. 


In  order  to  investigate  the  coalescence  question  further,  two  additional  studies  have 
been  carried  out.  First,  we  have  rezoned  the  calculations  to  allow  them  to  run  longer. 
Figure  14  shows  the  run  displayed  in  Figure  4,  run  to  longer  times  (w mt  =  16)  in  a  bigger 
mesh  ( Lx,Ly  — *  40 c/we).  Three  points  should  be  noted.  First,  after  the  initial  coalescence 
by  about  a  factor  of  two,  there  is  no  further  merging  of  the  modes.  Second,  while  the 
central  core  of  the  plasma  collapses  (along  with  the  magnetic  cavity),  the  outer  parts  of 
the  plasma  continue  to  expand,  as  was  noted  in  the  previous  section.  Third,  the  tips  of 
the  larger  structures  appear  to  bifurcate.  This  may  be  another  nonlinear  effect,  or  more 
likely,  is  the  reappearance  of  the  original,  shorter  wavelength  modes  [Winske,  1987], 

We  have  also  carried  out  a  detailed  study  of  the  mode  structure  for  a  selected  number 
of  runs.  To  show  the  propagation  of  the  waves  in  more  detail,  Figure  15  plots  the  density 
perturbations  as  a  function  of  6  at  the  front  edge  of  the  expanding  plasma  shown  in  Figure  4 
at  various  times.  The  average  value  of  n,  has  been  subtracted  out  and  the  amplitudes  have 
been  normalized  to  unity.  At  early  times,  one  sees  the  appearance  of  the  short  wavelength 
modes  which  move  slowly,  but  perceptively  in  the  positive  9  direction  (the  direction  of 
the  electron  drift).  Although  the  phase  velocities  of  these  waves  are  small,  the  fact  that 
k  is  large  (mode  ~  17  dominates)  allows  wr  to  be  appreciable.  In  fact,  from  measuring 
the  phase  speed  of  several  of  the  perturbations,  one  finds  u>r  ~  I.Swlh,  in  agreement  with 
linear  theory.  At  later  times,  these  shorter  wavelength  modes  merge  into  longer  wavelength 
structures  with  much  lower  phase  velocities  and  ur  ~  0. 

To  show  the  properties  of  the  waves  in  more  detail,  Figure  16  displays  the  Fourier 
spectra  of  the  density  perturbations  |n*|2,  k=l,2,...40  at  various  times.  The  spectra  are 
determined  from  Fourier  analyzing  the  density  fluctuations  every  0.25w~j ^  and  then  av¬ 
eraging  the  spectrum  over  several  intervals.  At  early  times  [uLjjt  =  2  -  3)  we  see  the 
emergence  of  short  wavelengths  (modes  m=13-17),  which  shift  to  slightly  higher  mode 
number  in  the  next  frames.  At  later  times,  longer  wavelength  modes  (m=7-8)  then  grow 
up  and  dominate.  This  can  be  seen  more  clearly  in  the  next  picture,  Figure  17,  where  the 
Fourier  modes  are  plotted  versus  time.  Again,  the  growth  of  shorter  wavelength  modes 
is  evident  for  u>nft  <  5,  with  the  longer  modes  dominating  later.  Note  that  there  does 
not  seem  to  be  a  cascade  through  a  progression  of  modes  down  to  m=7-8,  as  modes  9-12 
do  not  ever  grow  appreciably.  For  comparison  with  linear  theory,  the  solid  curve  denoted 
the  maximum  linear  growth  rate  (which  occurs  at  short  wavelength),  while  the  dashed 
curve  is  the  predicted  linear  growth  rate  for  modes  7-8.  The  results  suggest  that  the  late 
time  growth  of  modes  7-8  is  a  nonlinear  consequence,  not  a  slower  linear  growth  to  larger 
amplitude. 

Finally,  another  possible  nonlinear  effect  should  be  mentioned.  Sgro  et  al.  [1988] 
observed  no  coalescence  in  their  hybrid  simulations  with  m,/me  =  100,  but  a  coalescence 
by  about  a  factor  of  two  with  m,/me  =  1836.  While  the  results  are  consistent  with  those 
of  Figure  13,  Sgro  et  al.  [1988]  have  suggested  that  the  coalescence  scales  as  (m,/me)a 
with  a  ~  1/4.  (This  is  also  consistent  with  Figure  5.)  This  scaling  is  consistent  with  their 
simulations  and  when  scaled  to  the  mass  of  barium  ions  ( mBa/mp  =  137)  also  gives  the 
correct  wavelength  for  the  structures  observed  in  the  AMPTE  releases. 
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Figure  14.  Ion  density  contours  at  various  times  for  the  same  run  as  shown  in  Figure  4 
extended  to  later  time  by  expanding  the  simulation  box. 


Figure  15.  Density  perturbations  at  the  surface  of  the  expanding  plasma  versus  the 
imuthal  angle  0  at  various  times  for  the  run  shown  in  Figure  4. 


modes  7-8 


modes  13  —  14 


Figure  17.  Time  histories  of  the  Fourier  modes  shown  in  the  previous  figure.  Solid  lines 
correspond  to  the  maximum  linear  growth  rate;  dashed  line  is  the  linear  growth  rate 
for  modes  7-8. 


F.  Asymmetric  Expansions 

While  the  perfectly  symmetric  expansions  are  convenient  for  theory  and  numerical 
simulations,  one  expects  asymmetries  to  occur  in  laser  initiated  plasma  expansions  as  well 
as  in  VHANEs.  The  question  naturally  arises  whether  nonuniformities  in  the  expansion 
give  rise  to  new  effects.  Figure  18  shows  the  results  of  three  2-D  electromagnetic  particle 
simulations,  carried  out  to  late  times  by  expanding  the  grid,  at  several  points  in  time. 
The  top  expansion  is  uniform:  the  y-componeni  of  the  initial  expansion  velocity  Vy  =  V*, 
the  initial  x-component.  (The  particles  expand  radially  with  velocity  Vr  =  Vz  cos  Q  + 
V„sin0.)  The  middle  panels  correspond  to  Vy  =  0.5V'*,  the  bottom  panels  to  Vy  =  0.25VZ. 
Except  for  the  anticipated  fact  that  the  plasma  is  more  elongated  in  the  asymmetric 
expansions,  there  is  little  difference  in  the  mode  structure  at  early  times  and  the  amount 
of  nonlinear  coalescence  to  longer  wavelengths  is  about  the  same  in  each  case.  Of  course, 
the  asymmetries  (i.e.,  jets)  seen  in  the  laser  experiment  are  a  much  later  manifestation 
and  could  likely  imply  other  effect's,  but  at  early  times,  one  does  not  expect  the  short 
wavelength  physics  to  be  affected  very  much  by  asymmetries. 


G.  3-D  Simulations 

We  have  begun  a  study  of  3-D  effects  associated  with  expanding  plasmas  using  a 
newly  developed  extension  of  the  explicit  electromagnetic  particle  code  ISIS.  Here  we 
briefly  discuss  two  aspects  of  the  problem:  (l)  a  comparison  of  the  structuring  instability 
in  2-D  and  3-D,  and  (2)  a  new  effect,  streaming  of  energetic  electrons  from  the  cloud  along 
the  ambient  magnetic  field  [Barnes,  Jones,  and  Winske,  in  preparation]. 

The  2-D  simulations  are  similar  to  those  described  in  Section  2.  The  plasma  consists 
initially  of  a  thin  cylinder  of  radius  r0  =  2c/we.  The  electrons  and  ions  both  expand 
radially  with  velocity  Vr  =  0.0136c,  with  an  electron  (and  ion)  thermal  speed  of  0.001c, 
into  a  magnetic  field  B  =  Baz  of  strength  ne/we  =  0.8.  10000  particles  on  a  50  x  50  grid 
were  used  in  the  calculation.  The  left  panels  of  Figure  19  show  the  ions  in  the  x-y  plane 
at  two  times,  u >lhI  =  4.8  and  9.6.  As  with  the  simulations  described  earlier,  one  sees  a 
thin  plasma  shell  forms  and  large  flute  modes  grow  on  the  surface. 

The  3-D  simulation  starts  with  a  sphere  of  identical  radius,  with  the  particles  expand¬ 
ing  with  velocity  Vr  =  0.02.  The  radial  velocity  has  been  adjusted  so  that  Rb  is  the  same 
in  both  cases  [Eq.  (7)].  In  this  case  400000  particles  are  used  on  a  50  x  50  x  100  grid. 
The  right  panels  show  the  projection  of  the  ions  onto  the  x-y  plane  at  the  same  two  times. 
While  the  overall  effect  of  the  instability  is  qualitatively  the  same,  several  differences  be¬ 
tween  2  and  3-D  are  evident  from  the  figure.  First,  the  size  of  the  flutes  in  3-D  are  larger, 
suggesting  a  faster  growing  instability.  Second,  the  number  of  flutes  in  the  3-D  case  is 
slightly  less.  Both  of  these  features  are  consistent  with  the  larger  effective  deceleration 
in  3-D,  and  the  fact  thot  p,/Rb  is  slightly  larger.  Third,  the  3-D  case  seems  to  lack  a 
central  cavity;  this,  however,  is  a  visual  effect  resulting  from  projecting  all  the  ions  onto 
one  plane.  The  fact  that  the  flutes  are  visible  in  spite  of  the  projection  also  indicates  that 
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Figure  19.  Electromagnetic  particle  simulations  showing  ions  in  x-y  plane:  (left  panels) 
2-D  run;  (right  panels)  3-D  run. 


the  structures  are  field  aligned.  Finally,  there  is  a  marked  asymmetry  in  the  3-D  case. 
Such  asymmmetries  sometimes  also  occur  in  2-D,  and  it  is  not  thought  to  be  a  peculiarity 
of  3-D. 


However,  a  new  effect  is  observed  in  3-D,  as  shown  in  Figure  20.  Plotted  in  the  figure 
are  the  ions  and  electrons  projected  onto  the  y-z  plane  at  ujujt  =  4.8.  Again,  several 
interesting  features  are  visible.  First,  the  shape  of  the  ion  cloud  and  the  electron  cloud 
>  are  different.  The  ion  cloud  is  fatter  (in  y),  which  results  from  the  fact  that  the  ions 

expand  slightly  ahead  of  the  electrons  (recall  Figure  1).  This  is  referred  to  as  the  “ion 
charge  layer”  by  Galvez  et  al.  [1988].  Second,  the  cloud  is  not  spherical.  As  the  expansion 
stops  in  the  radial  direction,  the  motion  of  the  plasma  continues  along  z.  At  later  times, 
the  cloud  becomes  even  more  cigar  shaped,  as  was  observed  to  occur  in  the  AMPTE 
'  barium  releases  [Bernhardt  et  al.,  1987].  Third,  and  most  significant,  a  small  number  of 

electrons  have  streamed  out  of  the  cloud  along  the  magnetic  field.  This  effect  is  not  seen 
in  two-dimensional  simulations  in  which  B  is  in  the  plane  of  the  calculation. 

Two  tentative  explanations  for  this  phenomenon  are  as  follows.  It  is  observed  even  in 
I  the  2-D  simulations  where  the  magnetic  field  is  perpendicular  to  the  plane  of  the  calcu¬ 

lations,  that  in  addition  to  the  lower  hybrid  drift  instability  that  produces  the  structure 
on  the  surface,  another  instability  occurs  very  early  in  the  run.  This  instability  is  only 
seen  when  the  electrons  are  very  cold,  has  short  wavelength,  and  serves  to  heat  the  elec¬ 
trons  to  modest  temperatures  very  rapidly.  (Note  from  the  discussion  of  Sec.  3C  that  the 
|  slower  growing  lower  hybrid  drift  instability  is  hardly  affected  by  the  magnitude  of  Tt.) 

It  is  thought  that  the  instability  is  some  form  of  the  Kelvin-Helmholtz  instability,  similar 
to  that  observed  by  Sydora  et  al.  [1983]  in  their  electrostatic  simulations.  The  effect  of 
this  instability  in  3-D  is  to  allow  some  of  the  heated  electrons  to  escape  as  a  burst  along 
the  magnetic  field.  And  because  the  instability  does  not  occur  in  2-D  simulations  where 
|  B  is  in  the  plane  of  the  calculation,  the  electron  streaming  is  not  observed  in  that  case. 

A  second  possibility  is  that  the  lower  hybrid  waves  heat  the  electrons.  Even  though  the 
structures  appear  to  be  flute  modes  (k\\  =  0),  it  takes  only  a  small  component  of  k  along 
B  to  accelerate  the  electrons  [e.g.,  Tanaka  and  Papadopoulos,  1983].  More  study  is  needed 
to  distinguish  between  these  two  mechanisms.  Obviously,  because  the  plasma  will  try  to 
I  maintain  quasineutrality,  the  amount  of  electrons  which  can  escape  is  small.  Nevertheless, 

it  can  be  an  important  effect  and  one  which  we  will  investigate  more  completely  in  the 
future. 


I  H.  Weak  Instability  Limit 

We  have  seen  in  Sec.  3D  that  a  strong  instability  occurs  when  Pi/Re  >  1.  However, 
the  instability  can  still  be  excited  for  Pi/Ra  <  1.  When  g  =  0,  the  properties  of  the  lower 
hybrid  drift  instability  are  well  known  in  this  weak  drift  limit.  The  instability  persists 
I  down  to 

Ve/vi  =  4  (me/m<)1/4  (19) 


16 


40.0 

ions 

40.0 

electrons 

- t - t  r 

I  i 

y  . 

30.0 

— 

30.0 

20.0 

■  #  - 

x20.0 

-  #  - 

10.0 

— 

10.0 

i  •* 

0.0 

_ 1 _ i _ l . 

0.0 

.  i  _l 

-10.0  0.0  10.0  20.0  30.0  -10.0  0.0  10.0  20.0  30.0 

v  y 


Figure  20.  Results  of  3-D  simulations  showing  particles  projected  into  y-z  plane:  ions  (left 
panel),  electrons  (right  panel). 


Below  this  value  the  unmagnetized  ion  approximation  breaks  down  (w  <  ft,)  and  ion 
cyclotron  effects  stabilize  the  instability  [Freidberg  and  Gerwin,  1977], 

When  g  ^  0,  but  /?,  1,  the  instability  can  exist  at  even  very  weak  drifts.  Figure 

21  shows  the  results  of  linear  theory  with  cold  electrons,  unmagnetized  warm  ions  with 
enc/ui  =  2  and  m,/me  =  1836.  Plotted  are  the  growth  rates  maximized  over  k  as  a 
function  of  /?,  for  several  values  of  Vg/vA.  Generally,  7  decreases  with  both  Vg/vA  and 
Pi.  But  the  instability  still  occurs  as  Pi  — *  0  for  Vg  ±  0  and  close  to  pi  =  0  for  Vg  = 
0.  (The  bottom  curve  ends  where  uiT  <  ft,,  indicating  the  breakdown  of  the  theory.) 
Thus,  even  at  /?,  =  0.001,  for  Vg/vA  =  0.002,  implying  pi/R-B  =  0.0014  (for  a=1.5), 
one  finds  7  ~  0.08ft,.  And  the  number  of  growth  times  to  the  end  of  the  expansion  is 
7<  =  (7/fti)(ftjf)  (7/fti)(l-5f2fl/pv)  ~  80.  However,  as  we  shcrw  later  in  this  report, 
even  though  the  growth  rate  is  large  and  the  instability  has  plenty  of  time  to  grow,  its 
effect  on  diffusion  of  the  magnetic  field  and  heating  of  the  plasma  is  very  weak.  While  the 
instability  is  relatively  unimportant  in  this  limit,  one  must  also  keep  in  mind  that  other 
long  wavelength,  kpi  <c  1,  slower  growing  modes  probably  also  occur. 

I.  Collisional  Effects 

The  effect  of  collisions  with  the  background  plasma  is  important  in  the  laser  experi¬ 
ments  and  affects  HANEs  at  lower  altitudes.  The  study  of  collisional  processes  with  respect 
to  the  usual  lower  hybrid  drift  instability  has  been  investigated  in  some  detail  for  iono¬ 
spheric  applications  [Sperling  and  Goldman,  1980;  Huba  and  Ossakow,  1981;  Gary  et  al., 
1983],  Here  we  carry  out  a  simple  analysis  when  Vg  ^  0,  in  fact  Vg  >  Vn,  and  show  that 
such  effects  persist  and  can  be  important. 

The  inclusion  in  the  linear  theory  of  collisions  with  a  background  plasma  and  neutrals 
is  straightforward.  One  replaces  w  in  the  electron  terms  with  w  +  »i/e  and  with  w  +  t’17  in  the 
ion  terms.  Here  17.  =  uen  +  ueek2p2,  where  vtn  =  electron  neutral  collision  frequency  and 
i/ee  =  the  electron-electron  collision  frequency;  similar  expressions  occur  for  17  [Sperling 
and  Goldman,  1980].  In  addition,  the  presence  of  17  increases  the  drag  of  the  ions;  hence, 
the  magnitude  of  the  radial  electric  field  and  thus  Vg: 

cEr  —  T,'Cn  —  Tflil/iVi  (20) 

which  implies  an  increased  deceleration  g  =  -dVi/dt  +  17K,.  Figure  22  displays  growth 
rates  versus  k  for  the  same  parameters  as  Figure  3,  on  the  same  scale  for  easy  comparison. 
Here  we  have  assumed  Te  —  0  and  included  only  the  electron-neutral  and  iou-neutral 
collisions.  The  addition  of  ve  ^  0  reduces  the  maximum  growth  rate  but  does  not  shift 
the  wavelength  corresponding  to  the  most  unstable  mode  (as  would  occur  with  Tt  ±  0). 
Even  with  17  =  100ft,  ~  27  the  maximum  linear  growth  rate  is  not  reduced  very  much. 
However,  the  addition  of  ion  collisions  (here  17  =  20ft j)  has  a  marked  effect  on  the  growth 
of  the  instability  [Gary  et  al.,  1983]. 


iguxe  21.  Results  of  linear  theory  with  cold  electrons  showing  growth  rates  maximized 
over  wavenumber  versus  /?,  for  various  values  of  Vg/vA.  The  bottom  curve  cuts  off 
when  ur  <  indicating  the  breakdown  of  the  assumption  of  unmagnetized  ions. 


Figure  22.  Results  of  linear  theory  showing  growth  rates  as  a  function  of  wavenumber 
for  same  parameters  as  Figure  3,  except  that  vt  =  1000,  and  vt  =  0  (top  curve), 
i*i  —  200 *  (bottom  curve). 


The  fact  that  the  instability  persists  for  very  large  values  of  ut  is  well  known  [Gary 
et  al.,  1983].  Figure  23  shows  growth  rates  maximized  over  wavenumber  versus  ue/Qt 
for  various  values  of  Vg/v*  (with  Vn/vx  —  1,  /?,  =  0.2,  Tt  —  0,  17  =  0).  With  Vg  =  0 
(the  usual  lower  hybrid  drift  instability)  the  inclusion  of  collisions  strongly  reduces  the 
maximum  growth  rate,  but  does  not  change  the  wavenumber  corresponding  to  maximum 
growth  very  much.  An  extension  of  the  initial  fall  off  of  the  curve  with  ut  would  imply 
stability  when  ue  ~  7 (i/e  =  0).  However,  another  branch  of  the  instability  (the  collisional 
density  drift  instability  of  Gary  et  al.  [1983])  takes  over  and  gives  significant  growth  for 
ut  »  A,.  At  Vg  /  0,  these  two  distinct  modes  merge  into  one,  with  7  falling  off  slowly 
with  vt. 

The  effect  of  collisions  on  the  instability  has  been  tested  by  simulations  using  VENUS, 
a  2-D  electromagnetic  particle  code  which  allows  collisional  interactions  [Cranfill  et  al., 
1986].  In  this  preliminary  study  Rutherford  scattering  between  the  expanding  plasma 
particles,  rather  than  Krook  collisions  with  a  background  gas  that  was  assumed  in  the  linear 
theory  was  used.  Nevertheless,  the  effect  that  the  instability  may  persist  for  relatively 
high  rates  of  collisions  (relative  to  the  linear  growth  rate)  appears  to  occur.  Figure  24 
shows  snapshots  of  four  runs  at  one  instant  of  time.  For  each  run  the  ions  in  x-y  space, 
contours  of  ion  density,  and  magnetic  field  contours  are  shown.  From  top  to  bottom  the 
collision  frequency  was  varied  from  ue  =  0, 0.1, 10, 100a >lh  (t'e/17  ^  30).  For  large  collision 
frequencies  the  plasma  does  not  expand  very  much,  as  the  directed  motion  is  quickly 
thermalized  with  only  a  few  energetic  particles  escaping.  Nevertheless,  an  instability  at 
roughly  constant  wavelength  appears  in  each  case.  Obviously,  much  more  work  using  a 
collision  model  consistent  with  the  theory  needs  to  be  carried  out,  but  the  initial  results 
are  encouraging. 


4.  Discussion 

This  final  section  consists  of  three  parts:  (A)  a  discussion  of  issues  regarding  the 
theoretical  work  and  its  relation  to  the  laser  experiments,  (B)  the  consequences  of  the 
results  for  very  high  altitude  nuclear  explosions,  and  (C)  a  short  summary. 

A.  Laser  Experiments 

The  results  presented  in  this  report  raise  a  number  of  questions  regarding  the  find¬ 
ings  of  the  several  laser  experiments  in  which  structures  on  expanding  plasmas  have  been 
observed.  We  subdivide  the  subsequent  discussion  into  four  areas:  (1)  properties  of  the 
observed  waves,  (2)  pi/ Rb  effects,  (3)  late  time  behavior,  and  (4)  collisional  effects. 
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Figure  24.  Results  of  2-D  electromagnetic  particle  simulations  with  Rutherford  collisions: 
ions  in  real  space,  ion  density  contours,  and  magnetic  field  contours  at  one  time  for 
four  runs  with  different  values  of  the  collision  frequency. 


1.  Observed  Waves 


The  key  question  is  what  is  the  relation  of  the  waves  observed  in  the  simulations  to 
those  seen  in  the  experiments.  It  has  been  argued  that  the  lower  hybrid  drift  instability 
cannot  be  responsible  for  the  structures  seen  on  the  surface  because  the  theoretical  wave¬ 
lengths  are  much  too  small  compared  to  those  measured  in  the  NRL  experiments  [Ripin  et 
al.,  1987]  and  in  the  AMPTE  magnetotail  releases  [Bernhardt  et  al.,  1987].  Such  estimates 
are  usually  based  on  the  expression  kpe  ~  1  (Te  7V)  and  assuming  a  very  cold  electron 
temperature.  However,  we  have  seen  that  important  linear  modifications  including  g  ^  0 
(Sec.  2),  correct  inclusion  of  3-D  dynamics  (Sec.  3A),  finite  0e  (Sec.  3B)  and  nonlinear 
effects  (Sec.  3D)  all  operate  to  push  the  theoretical  value  of  k  to  smaller  values.  There 
are  also  other  effects,  e.g.,  nonlocal  theory  with  electromagnetic  corrections  and  realistic 
radial  profiles,  which  would  also  modify  the  wavelength  of  the  instability.  A  related  dif¬ 
ficulty  with  comparing  with  experiments  is  that  often  the  parameters,  such  as  the  local 
value  of  the  deceleration  and  the  ion  temperature,  are  not  known  with  much  certainty.  In 
addition,  other  effects,  such  as  the  changing  charge  state  of  the  aluminum  plasma  in  the 
NRL  experiment,  that  are  not  included  in  the  theory  as  yet,  may  also  be  important. 

A  second  issue  involving  the  wavelength  of  the  surface  flutes  is  how  they  scale  with 
the  applied  magnetic  field,  Bc.  While  the  Japanese  experiments  show  k  ~  B°-8,  the  NRL 
experiments  do  not  show  any  consistent  variation  of  k  with  B„.  One  possible  explanation 
is  suggested  by  the  simulations  and  was  shown  in  Figure  5.  Measurements  of  k  during 
the  linear  growth  phase  indicated  k  ~  B°-5,  consistent  with  linear  theory  when  k  was 
normalized  to  Rb-  However,  during  the  nonlinear  phase  when  the  modes  coalesced,  k  was 
roughly  independent  of  B0.  It  may  be  that  in  the  NRL  experiments  (at  larger  Pi/Rb) 
that  one  is  observing  the  structures  in  their  nonlinear  state,  rather  then  while  they  are 
growing  exponentially  according  to  linear  theory,  which  may  be  more  characteristic  of  the 
Okada  et  al.  [1981]  observations.  This  question  may  be  resolved  when  the  NRL  experiment 
acquires  the  capability  to  make  multiple  observations  per  shot.  In  a  similar  vein,  it  would 
be  useful  to  examine  the  variation  of  the  wavelengths  of  the  flutes  with  the  ion  to  electron 
mass  ratio,  by  changing  the  laser  target  material.  Figure  5  suggests  a  weak  (m./me)1^4 
dependence  in  both  the  linear  and  nonlinear  regimes,  as  do  the  simulations  of  Sgro  et  al. 
[1988]. 

A  related  issue  involves  the  real  frequency  of  the  waves.  Both  the  laser  experiments 
and  the  simulations  show  flutes  that  appear  to  be  stationary,  implying  ur  —  0.  However, 
the  most  unstable  modes  according  to  linear  theory  are  characterized  by  7  ~  ujt  ~  ulh- 
The  simulations  do  show,  however,  that  by  the  time  the  flutes  are  visible  in  the  particles 
(which  would  correspond  to  optical  observations)  that  the  waves  are  already  rather  non¬ 
linear  (e6<f>/Ti  >  1).  At  early  times  the  frequency  is  large  (<*v  ~  wlh),  but  the  phase 
velocities  are  small  because  k  is  large.  At  later  time  the  nonlinearities  produce  phase  shifts 
of  the  frequency  so  that  wr  — *  0.  Indeed  such  phase  shifts  would  allow  the  mode  coupling 
process  to  proceed  more  easily  by  allowing  the  k’s  and  u/’s  to  satisfy  the  matching  condi¬ 
tions,  k0  =  ki  +  *2i  w0  =  wi  +W2-  Such  nonlinear  frequency  shifts  are  well  known  to  occur 
for  the  lower  hybrid  drift  instability  with  g  =  0  [Winske  and  Liewer,  1978]. 
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All  the  laser  experiments  [Ripin  et  al.,  1987;  Okada  et  al.,  1981;  Zakharov  et  al.,  1986] 
also  detected  high  frequency  plasma  noise  in  the  lower  hybrid  frequency  range.  Similar 
noise  was  measured  by  the  AMPTE/IRM  satellite  (Gurnett  et  al.,  1986]  at  the  edge  of 
the  cavity  in  both  the  electric  and  magnetic  receivers.  Because  the  lower  hybrid  drift 
instability  is  essentially  an  electrostatic  mode,  but  also  has  an  electromagnetic  component, 
it  is  the  likely  source  of  the  noise.  However,  since  the  instability  is  not  sharply  peaked  in 
frequency,  one  expects  a  broad,  rather  featureless  spectrum  is  generated,  consistent  with 
the  measurements.  Even  though  we  have  previously  argued  that  the  dominant  modes  may 
nonlinearly  shift  their  frequency  to  smaller  values,  a  range  of  low  amplitude  waves  will 
still  be  continuously  generated  and  provide  a  broad  frequency  spectrum.  Although  the 
spectrum  is  broad,  the  fact  that  the  frequency  centers  on  ulh  suggests  looking  for  shifts 
in  the  bandwidth  with  the  applied  magnetic  field  and  the  mass  of  the  ions. 

A  final  issue  concerns  the  characteristics  of  the  observed  structures  in  asymmetric 
expansions.  Our  simulations,  which  describe  the  early  time  behavior,  suggest  that  asym¬ 
metries  in  the  expansion  do  not  affect  the  properties  of  the  short  wavelength  modes  nor 
their  coalescence  to  larger  structures  appreciably.  The  principal  effect  would  come  from 
a  modified  effective  deceleration  (g)  which  an  asymmetry  could  produce.  A  study  of  the 
wavelengths  of  the  observed  flutes  as  the  asymmetry  is  varied  (by  changing  the  angular 
spread  of  the  laser)  could  shed  future  light  on  this  question. 


2.  Px/Rb  Effects 

Various  experiments  have  shown,  and  the  simulations  here  have  confirmed,  that  the 
single  most  important  parameter  is  p,/Rb-  When  Px/Rb  >  1,  as  in  the  NRL  experiments, 
large  flutes  are  observed  that  display  interesting  nonlinear  dynamics  and  the  size  of  the 
magnetic  cavity  (Rm)  is  less  than  RB  [Zakharov  et  al.,  1986].  When  Pi/Rs  <  1,  as  in  the 
AMPTE  releases,  the  Japanese  laser  experiments,  and  theta  pinches,  the  flutes  are  smaller 
and  exhibit  less  dramatic  behavior,  while  Rm  a  Rb.  When  Px/Rb  <■  1,  as  in  the  hot 
diamagnetic  cavities  upstream  of  the  Earth’s  bow  shock  (and  probably  also  a  VHANE), 
only  a  very  weak  instability  is  detected.  Such  results  are  also  consistent  with  theory  in 
that  Px/Rb  ~  Vfl/tu-  For  moderately  cold  plasmas  the  contribution  to  the  cross-field  drift 
due  to  Vg  is  greater  than  that  due  to  the  usual  diamagnetic  drift  Vn,  so  that  Vg  is  a  direct 
measure  of  the  free  energy  of  the  instability. 

A  most  useful  set  of  laser  experiments  would  be  to  vary  Px/Rb  systematically.  In 
2-D,  as  in  the  simulations,  both  Rb  and  pi  vary  as  B~l,  so  Px/Rb  is  independent  of 
B0 ;  however,  in  3-D  RB  ~  Bo*^3>  so  pi/Rs  ~  Bo1/3.  Other  ways  to  change  Px/Rb 
would  be  to  either  decrease  the  laser  energy  so  that  the  streaming  velocity  VD  is  reduced, 
Px/Rb  ~  Vp  l/Z,  or  change  the  target  material  to  alter  the  ion  mass.  One  would  like  to 
quantify  the  wavelengths  and  the  radial  extent  of  the  flutes  (relative  to  Rb).  A  more 
difficult  measurement  would  be  the  size  of  the  magnetic  cavity  compared  to  Rb,  again  for 
a  range  of  Px/Rb  [e.g.,  Zakharov  et  al.,  1986]. 
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3.  Nonlinear  Effects 


The  surface  modes  generated  in  the  pxf  Rb  >  1  regime  exhibit  a  number  of  interesting, 
probably  nonlinear  effects.  These  include  the  freestreaming  of  the  flute  tips,  bending  of 
the  flutes  in  the  magnetic  field,  and  bifurcation  of  the  ends  of  the  structures  [Ripin  et  al., 
1987]. 

The  simulations  carried  out  in  this  regime  show  a  number  of  similar  features  [Winske, 
1987].  Continued  outward  streaming  of  the  flutes  is  seen,  while  the  inner  portion  of  the 
plasma  and  the  magnetic  cavity  recollapse  [e.g.,  Figures  4  and  12].  Bending  of  the  flutes  is 
often  observed,  although  the  sense  of  the  bending  is  usually  opposite  to  the  experiments. 
The  flutes  in  the  simulations  bend  clockwise,  following  the  ion  sense  of  rotation,  rather  the 
electron  sense.  An  exception  are  the  hybrid  simulations  of  Sgro  et  al.  [1988],  where  the 
flutes  bend  in  the  electron  sense.  In  this  case,  however,  the  current  seems  to  be  carried 
mainly  by  the  ions,  which  have  a  large  diamagnetic  drift. 

Furthermore,  the  simulations  clearly  show  mode  coalescence  and  later  the  reappear¬ 
ance  of  shorter  wavelength  modes  [e.g.,  Figure  14].  It  appears  from  am  amalysis  of  the 
simulations  that  these  short  wavelength  structures  are  remnants  of  the  original  linear 
modes,  rather  than  bifurcation  of  the  larger  structures.  It  would  be  interesting  in  both 
the  simulations  and  the  laser  experiments  to  check  if  these  larger  k  modes  exhibit  some 
dependence  on  fl0,  as  do  the  waves  in  the  linear  regime. 


4.  Collisional  Effects 

A  final  effect  that  we  have  studied  involves  collisions  with  the  background  plasma. 
Both  linear  theory  and  the  particle  simulations  suggest  that  the  instability  can  persist  for 
very  large  values  of  the  electron-neutral  collision  frequency  at  roughly  the  same  wavelength. 
On  the  other  hand,  ion  collisions  serve  to  both  increase  the  effective  deceleration  (which 
increases  the  growth  rate  slightly)  and  more  importantly,  to  damp  the  waves. 

In  the  laser  experiments  the  collision  frequency  can  be  varied  by  changing  the  fill 
pressure  of  the  background  gas.  The  experiments  at  NRL  suggest  that  as  the  fill  pressure 
is  raised,  the  instability  is  less  easy  to  see  optically,  retains  about  the  same  wavelength,  and 
is  eventually  quenched  if  the  background  pressure  is  too  high.  Overall,  these  qualitative 
results  are  consistent  with  the  simple  theoretical  analysis  that  has  been  done  so  far.  A 
more  careful  study  of  this  effect  experimentally  both  as  a  function  of  fill  pressure  and  if 
possible  Te/Ti,  to  change  Vt/ui,  would  be  useful. 

B.  Consequences  for  VHANE 

We  next  turn  to  a  very  different  parameter  regime,  that  of  a  very  high  altitude  nu¬ 
clear  explosion  (>  1000  km).  Assuming  a  weapon  yield  of  about  one  megaton  with  about 
a  quarter  of  the  energy  going  into  kinetic  energy  of  the  debris  ions,  yields  about  N  ~  1028 
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ions.  Assuming  an  expansion  velocity  of  2  x  103km/sec  and  a  simple  pressure  balance  rela- 
■  tionship  (7)  with  B0  =  1/3 G  gives  Rb  ^  800km.  Taking  A  =  27,  Z  =  13  for  an  “average” 

ion,  one  finds  p,-  ~  1.3  km.  Thus,  Pi/RB  =  1.6  xlO-"1  and  V9/va  —  1.5 Pi/Ra  =  2.4xl0~3, 
in  contrast  to  the  laser  experiments  and  simulations  where  Pi/ Rb  ~  1  Furthermore,  tak¬ 
ing  the  plasma  temperature  ~  10  keV  yields  Pi  =  0.002  at  maximum  expansion,  which  for 
tnc/ui  ~  2  gives  Vn(vx  =  0.002.  Thus,  although  one  has  Vb/va  C  1  and  Vnjvx  <■  1,  still 
|  Vn  ~  V„,  as  in  the  simulations. 

Even  though  the  gradients  and  cross-field  drifts  are  very  small  under  these  conditions, 
the  lower  hybrid  drift  instability  is  still  excited.  Figure  18  shows  growth  rates  maximized 
over  wavenumber  versus  Pi  in  this  regime.  Under  conditions  expected  for  a  VHANE,  the 
growth  rate  is  still  significant  7/flj  ~  0.1.  As  shown  in  Sec.  3H,  the  number  of  e-foldings 
I  in  this  regime  is  large,  so  that  7 1  ~  100.  But  the  important  question  remains,  how  large 

are  the  anomalous  effects  associated  with  the  instability  in  this  case. 

Figure  18  shows  that  for  Pi  small,  7  ~  ~  (V^/t;,)2.  This  weak  drift  limit  for  the 

usual  (no  g)  lower  hybrid  drift  instability  has  been  investigated  in  much  detail.  Both  the 
1  linear  properties  of  the  waves  [Davidson  and  Gladd,  1975]  and  their  nonlinear  consequences 

[Brackbill  et  al.,  1984]  are  well  understood.  Here  we  will  use  some  of  these  known  properties 
to  estimate  what  will  happen  in  a  VHANE.  Unfortunately,  the  PijRB  <-  1  regime  is  not 
readily  accessible  to  simulations,  as  the  instability  is  so  weak. 

|  The  anomalous  collision  frequency  for  the  lower  hybrid  drift  instability  is  given  by 

Liewer  and  Davidson  [1977]: 


k'an 


nmeV J/2 


(21) 


For  a  VHANE,  we  take  Ve  =  Vg  +  Vn  ~  2Vn;  (V^/v,)2  ~  4 (Vn/«i)2  ~  4/?,-  and  take 
€f  =the  nonlinear  fluctuation  level  of  the  waves  22  nmeVj/2  [Brackbill  et  al.,  1984].  At 
r  =  Rb  the  debris  density  is  n  ~  5  x  103  so  that  w2/D2  ~  1  and  uan  ~  0.0045wlh  ~  fl ,. 
Using  this  value  and  following  the  discussion  of  Sec.  3D  we  can  calculate  the  amount  of 
anomalous  diffusion  expected:  Ax2  ~  Dt  ~  $Vp{l2(pi/ RB)3t4  or  Ai /Rb  ~  10-4,  which 
is  negligible.  In  other  words,  the  plasma  expansion  is  essentially  classical  with  regard  to 
this  instability.  (Of  course,  other  types  of  slower  growing,  longer  wavelength  instabilities 
that  we  have  not  considered  here  can  still  be  operative.) 


In  similar  fashion  one  can  calculate  the  anomalous  heating  rate  for  the  ions,  Q'an  ~  uan 
[Liewer  and  Davidson,  1978] 


±d7\  &  meV 2 

Ti  dt  Von  3Ti 


2  x  10-9  Cti 


(22) 


So  negligible  ion  (as  well  as  electron)  heating  is  likely.  In  addition,  it  was  noted  in  Sec.  3G 
that  electron  acceleration  and  streaming  along  the  magnetic  field  may  occur  during  the 
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expansion.  However,  such  effects  require  very  cold  electrons  initially  and  perhaps  special 
initial  conditions,  as  well  as  significant  deceleration.  None  of  these  conditions  is  likely  to 
occur  early  in  a  VHANE.  Thermal  (i.e.,  ~  10  keV)  electrons  will  have  a  directed  velocity 
comparable  to  their  thermal  speed  and  for  px/Rb  <  1  the  deceleration  will  be  weak. 
Acceleration  of  electrons  by  lower  hybrid  waves  [e.g.,  Tanaka  and  Papadopoulos,  1983]  is 
also  possible,  and  seen  in  the  3-D  simulations,  but  again  is  not  likely  due  to  the  very  low 
level  of  fluctuations. 

Thus,  one  does  not  expect  the  lower  hybrid  drift  instability  to  be  very  important 
for  VHANES.  It  is  most  likely  that  the  cavity  expands  in  almost  a  classical  fashion  to  its 
maximum  radius,  with  a  low  level  of  short  wavelength  fluctuations  on  the  surface.  Because 
the  instability  is  expected  to  be  so  weak,  no  strong  nonlinear  effects  are  likely  to  occur: 
i.e.,  no  large  scale  flutes  are  expected  to  form  and  break  off.  And  because  the  diffusion  is 
very  small,  the  plasma  will  contract  radially  and  continue  to  expand  axially.  Of  course, 
the  large  size  of  the  cavity  and  how  it  interacts  with  the  Earth’s  diverging  magnetic  field 
has  yet  to  be  included  in  the  analysis.  But  such  processes  can  be  considered  in  the  absence 
of  strong  anomalous  processes  fairly  easily. 

That  such  large  scale  cavities  are  rather  benign  is  not  suprising.  Large  diamagnetic 
cavities  occur  upstream  of  the  Earth’s  bow  shock.  These  typically  are  larger  than  one 
Earth  radius  (Re  ~  6400km),  while  px  for  a  typical  solar  wind  proton  (VD  ~  500km, 
B  =  10~4G)  is  about  500  km  so  that  Px/Rb  <  0.1.  Although  one  is  limited  by  the  usual 
problem  of  separating  temporal  and  spatial  behavior  with  single  spacecraft  measurements, 
there  is  no  indication  that  large  scale  surface  structures  or  asymmetries  exist  [Thomsen 
et  al.,  1986],  Although  plasma  noise  is  detected  at  the  surface  of  the  cavity,  the  measured 
gradients  and  plasma  temperature  do  not  suggest  strong  anomalous  effects.  The  ground 
based  and  spacecraft  observations  of  the  cavity  produced  by  the  AMPTE  barium  releases 
in  the  magnetotail  in  which  pi /Re  ~  1  again  were  consistent  with  a  rather  nonviolent 
behavior.  [Bernhardt  et  al.,  1987].  The  observed  flutes  were  small,  did  not  seem  to  grow 
into  extended  fingers,  and  the  collapse  of  the  cavity  according  to  classical  behavior  again 
suggested  a  very  nondramatic  instability 

One  area  where  the  instability  may  be  relevant,  however,  is  to  lower  altitude  bursts, 
like  CHECKMATE,  at  very  early  times.  In  such  situations  collisional  effects  with  the  back¬ 
ground  can  be  important.  We  have  seen  that  the  lower  hybrid  drift  instability  can  persist 
in  the  presence  of  large  collision  frequencies  in  Sec.  31.  A  separate  report  in  preparation 
discusses  possible  application  of  the  lower  hybrid  drift  instability  in  this  more  collisional 
regime. 

C.  Summary 

To  sum  up  briefly,  in  this  report  we  have  discussed  a  number  of  issues  related  to  early 
time  structuring  in  very  high  altitude  nuclear  explosions  and  in  laser  experiments  that 
study  some  of  the  physics.  Our  emphasis  has  been  on  retaining  electron  effects  in  the 
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analysis  throughout.  The  inclusion  of  electron  physics  determines  a  scale  length  on  which 
the  fastest  growing  structures  develop.  We  have  shown  how  the  longer  wavelength  effects 
that  are  observed  can  arise  from  these  smaller  seed  structures.  We  have  reviewed  the  basic 
mechanism  of  the  instability,  the  results  of  a  linear  analysis  based  on  cold  electrons,  and 
particle  simulations  in  which  both  electron  and  ion  kinetics  are  included.  Furthermore,  we 
have  discussed  a  number  of  enhancements  to  the  model,  both  in  the  linear  theory  and  in 
the  interpretation  of  a  number  of  simulations  carried  out  over  a  range  of  parameters.  We 
have  also  related  the  results  to  the  interpretation  of  the  NRL  observations  and  suggested 
some  experimental  tests  of  the  theory.  Finally,  we  have  addressed  how  the  theory  scales 
to  the  VHANE  situation  and  indicated  areas  for  further  study  in  this  regime. 
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